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ABSTRACT 


Tn  Section  I the  Rayleigh  reflection  Coefficient  for  a rough 
surface  is  derived  in  a form  suitable  for  inclusion  in  the 
basic  integral  for  the  scattered  pressure.  A formula  for  the 
pressure  scattered  from  a rough  attenuating  surface  is  then 
derived  using  this  coefficient.  The  backscattered  intensity 
is  calculated  in  Section  TI  for  various  assumed  distributions 
of  surface  heights.  The  results  presented  here  are  restricted 
to  the  high  frequency  limit,  and  are  developed  using  both  the 
stationary  phase  and  the  modified  Fresnel  techniques.  Finally, 
a backscattering  result  for  a composite  surface  is  derived, 
and  is  shown  to  be  consistent  for  both  forward  and  backward 
scattering. 
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T . SCATTERED  PRESSURE  FROM  A RANDOMLY  ROUGH  ATTENUATING  SURFACE 

During  this  quarter  it  was  discovered  that  the  expressions  given 
in  the  previous  progress  report  for  the  Rayleigh  reflection  coefficient 
for  a rough  surface  were  inadequate.  The  results  given  previously 
were  only  approximate  due  to  an  error  in  the  method  in  which  Snell's 
law  was  applied  to  the  coefficient.  This  quarter  a more  complete 
and  accurate  treatment  of  this  problem  was  derived.  The  following 
derivation  differs  from  the  original  treatment  of  the  problem  not  only 
in  the  form  of  the  Rayleigh  coefficient  itself,  but  also  in  that  the 
coefficient  is  used  in  the  context  of  the  original  integral  for  the 
scattered  pressure. 

The  scattered  pressure  at  point  A is  given  by  the  integral 


where 

N is  a variable  point  on  the  surface, 
n is  the  outward  normal  to  the  surface,  and 

r^  is  defined  in  Dwg.  AS-68-1100.  The  exact  form  of  the  Kirchhoff 
boundary  condition  is 


where  p (N)  and  p are  the  values  of  the  scattered  and  incident 
pressures  on  the  surface  and  R is  the  Rayleigh  reflection  coefficient 
which  will  later  be  shown  to  be  a function  of  the  local  slope  of  the 
suiiace.  If  it  is  now  assumed  that  the  incident  pressure  is  of  the 
form 


ikr 
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it  is  then  obvious  that 
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p (N)  = RP  
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Following  the  development  given  in  Section  III  - D of  the  Final  Report 
under  Contract  1:00024-69-0-1275,  this  integral  becomes 
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where 


y - sinG,  + sin© 

1 r 


(7) 


und  is  a unit  vector  directed  along  line  R^.  Continuing  this 
derivation,  it  is  now  assumed  that 


A A 

e • e.  = cos 0 
y 1 r 


e • e = sin© 
z 1 r 
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so  that 
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xV5/y-y  • "l  = -Sln0r  + CyCOS0r 
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rhus  the  integral  for  the  scattered  pressure  is 
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sin0^+f ^cos0^)|  dxdy  . (10) 


Before  this  integral  can  be  further  simplified,  a suitable  expression 
for  the  Rayleigh  coefficient  must  be  derived.  Referring  to 
Dwg.  AS-7O-795,  the  Rayleigh  reflection  coefficient  is  given  by 


r(*1,«2) 


p„c  sin<J>  - p,c  sin<J>„ 

g-  g-  1 11 £ 

p2C2sin^]  + P1c1Einf>.. 


(11) 


where  ^ = 0^  + p and  t>,.  = 0,  + p 
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At  this  point  Snell's  law  must  be  applied  to  the  formula  by  making  the 
substitution 


sin<I>  = V 1 - cos"-$n  = -y  Vw'  - cos  $ , ( 1J ) 

d c.  JN  1 

where  N is  the  acoustic  index  of  refraction.  Applying  Eq.  (12)  and 
Eq.  (13)  to  Eq.  (ll)  gives 

p.  c r— 

p^c  sin(0  +p)  - — — Vnc  - cos~(0  +p) 

K(e1+P).^ 1 ^ — ■ dM 

P2c2sin(0^+P)  + - cosdq+P) 

By  expanding  sin^+p)  and  cos^+p)  and  using  various  trigonometric 
identities,  this  may  be  arranged  into  the  form 
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where 


q = tanp  , and 
A = pocosin0., 
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or  ease  of  application,  R(0^>7)  is  expanded  in  a Maclaurin  ceric r:: 


The  first  term  of  this  series  is  just  the  Rayleigh  reflection 
coefficient  for  a plane  surface. 

Since  7 is  the  tangent  of  p,  it  is  obvious  from  Dwg.  A3-70-797 
that  7 is  just  the  slope  of  the  locaL  tangent.  Because  the  formula- 
tion is  restricted  to  the  plane  of  incidence,  it  may  be  assumed  that 
7~ry.  Using  this  in  the  original  pressure  integral  gives 

r r / ih(  1-,  +K^  )\ 

P (a)  = II  p (- — ) e‘ik7  ( e , ,tj)(  -sine+nc  060  j]  is  ly 


At  this  point  £ and  r)  are  recognized  as  random  variables  and  the 

integral  for  p (A)  must  be  replaced  by  <p  (a)>.  If  it  is  assumed 
s s 

that  the  surface  height,  ?,  and  the  surface  slope,  tj  , are  independent, 
then  the  integral  for  <p^(A)>  takes  on  the  simple  form 

^.(e^X  _sin9^+T]cos0r 

(23) 

Since  it  has  been  assumed  that  tj  is  stationary  and  independent  of  t , 
it  may  be  shown  that 


<Ps(A)>  = ^R(01,T1)(l-rlcot0r 
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Since  ? is  also  stationary,  its  expectation  value  may  also  be  removed 
from  the  integral,  leaving  only  a form  which  is  immediately  recognized 
as  the  pressure  scattered  from  a plane  surface. 


<PS(A)>  = [<r(01,tj)>  - cot0r<nR(ei,Ti)>J^'lk"'r^i 
Assuming  a Gaussian  distribution  of  heights, 

_g 

<PS(A)>  = £<R(  Q-,  ,T))>  - coter<r,R(01,T)>J  e ' Pr 


(25) 


26 
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where 


g = ( k7o)", 

9 is  the  receiver  angle,  and 

p ^ is  the  pressure  scattered  from  a plane  surface. 

This  form  is  particularly  simple  to  calculate  due  to  the  form  of  the 
expansion  of  R(S  .t).  Rearranging  terms  of  the  series  gives 

|<K>  - cot0r<T]F>"j  = [R(l  - <ri>  cotf^  ) + R'(<r]>  - <Tj~>  cct$T  ) 


(<H  > - <r\J>  cot©  ) + 


(27) 


Thus  the  only  expectation  values  to  he  calculated  are  of  the  form 


00 

j ti4P(tj)cLti 

-tan# 


(28) 


where  the  lov:er  limit  of  -tan#^  represents  a simple  shadowing 
consideration.  For  the  Gaussian  case  these  integrals  are  quite  easy 
to  calculate  and  become 


where  s is  the  rrr.s  slope  of  the  unshadowed  surface. 
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For’  the  case  where  l is  an  odd  integer,  this  becomes 
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so  that 
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ihio.  may  be  easily  integrated  using  Gradshteyn  and  F.yzhiK  integral 
"O.  3-331-3  (p-  317,  ^th  edition)  and  is  found  to  be 
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where  P(a,x)  is  the  incomplete  gamma  function. 


For  the  case  where  i is  an  even  integer 


tanS. 
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These  values  may  he  easily  calculated  and  combined  with  the  values 
of  R,  4—  , —4  , etc.,  to  find  the  value  of  r<p>  - cot©  <nR>|  . 

^ 1=0  ^ ,„o  L ' J 

The  values  of  R,  R',  and  P"  will  be  complex  because  of  the  assumed 
attenuation  in  the  bottom,  so  the  entire  expression  for  the  pressure 
will  also  be  complex.  An  initial  numerical  calculation  was  made  using 
only  the  first  two  terms  of  the  expansion  of  R(0  ,tj)  which  gave  a 
reasonable  agreement  with  Mackenzie's  theoretical  curves,  which  were 
presented  in  the  previous  progress  report.  However,  the  lack  of  agree- 
ment in  certain  areas  of  the  curves  indicates  the  need  to  include 
higher  order  terms  in  the  calculation. 
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II.  BAC SCATTERING  THEORY  AND  EXPERIMENT 


Introduction 


The  analysis  of  experimental  backscattering  data  proves  much  more 
difficult  than  analysis  of  either  forward  or  specular  scattering  data. 
Many  investigators  have  assumed  a composite  (several  types  of  inde- 
pendent irregularities)  scattering  mechanism  in  order  to  obtain  the 
desired  fit  to  their  experimental  data.  However,  if  the  same  composite 
theory  (with  the  same  input  parameters)  is  applied  to  forward  scat- 
tering predictions,  satisfactory  agreements  are  not  obtained.  In 
addition,  many  composite  treatments  are  limited  to  use  of  the  nonana- 
lytic  exponential  correlation  function.  In  Section  D of  this  chapter 
a qualitative  description  of  a composite  theory  (two  types  of  roughness 
in  this  case)  is  given  which  does  not  have  these  limitations.  However, 
the  assumption  of  two  types  of  irregularities  requires  that  the 
surface  statistics  be  known  in  great  detail.  In  fact,  it  is  difficult 
to  see  how  the  composite  assumtion  is  going  to  help  the  backscattering 
prediction  problem  other  than  to  provide  an  "adjustable"  parameter. 


Section  B will  be  devoted  to  an  alternate  method  which  yiel 1. 
final  results  similar  to  those  obtained  by  the  composite  approach  by 
using  an  exponential  correlation  function.  The  method  used  will  be 
the  stationary  phase  technique  developed  in  the  Final  Report  under 
Contract  N00024-69-C -1275 . The  stationary  phase  equations  yield 
good  agreement  with  the  forward  and  spec- alar  scattering  data. 
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stationary  Phase  Backscatter 


The  "basic  formula  for  the  stationary  phase  method  is  taken  as 
Sq.  (49)  from  the  Final  Report  under  Contract  N00024-69-C -1275 . 
Equation  (49)  states  that  the  scattering  coefficient  a (not  the 
scattering  strength)  can  "be  written  as  follows: 
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where  the  notation  and  geometry  are  used  in  the  previously  referenced 
final  report.  Equation  (36)  is  next  integrated  for  the  case  of  a 
three-halves  distribution  of  heights. 


The  marginal  probability  density  function  for  the  three-halves 
distribution  of  heights  is  given  by 


P(’) 


and  the  bivariate  distribution  is  given  by 


p(V.')  = 


2*h?(l-C)V* 
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Tince  a variance  does  not  exist  for  these  distributions,  a convenient 
quantity  to  equate  to  the  experimental  rms  height  is  given  by  the 
interquartile  range.  The  interquartile  range  can  be  written  as  — h. 


The  characteristic  functions  associated  with  Eq.  (37)  and  Eq.  (38) 
can  be  stated  as 


/-ik7^  = krhK^(k7h) 

^-ik7(C-t')^  = e-  [2g(  1-C  )] 1 " 

2 2 p 

where  g = k 7 IT  and  h = 0.866  hg  ( experimental  rms  height). 

When  Eqs.  (39)  and  (4o)  are  used  in  Eq.  (36)  we  obtain  the 
following  scattering  coefficient: 
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The  asymptotic  form  of  Eq.  (4l)  can  be  stated  as 
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..here  a Gaussian  correlation  function  (with  correlation  length  L) 
has  been  assumed.  If  the  full  series  expansion  for  Eq.  (4o)  is 
applied  to  Eq.  (4l),  we  obtain  the  following  form: 
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where  the  Gaussian  correlation  function  again  has  been  used.  In 
order  to  compare  with  other  authors' work,  Eq.  (42)  can  also  be 
stated  ( for  backscattering)  as 


ahf  02gK^(g1/2)  + 
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where 

2 2 2 
C = L cos  «/ 2h  , 

F = l/cos®,  and 

represents  the  angle  of  incidence,  (i.e.,  $ + 6 = 90  deg). 

The  expression  in  braces  is  of  the  same  form  obtained  by  Fung  and 
Leovaris  [Ref.  l],  Beckmann  [Ref.  2],  and  Marcus  [Ref.  3].  However, 
these  authors  used  the  nonanalytic  exponential  correlation  function 


16 


J 


in  deriving  their  results.  As  recently  shown  by  Barrick  [iv.-f.  U:, 
the  use  of  the  exponential  correlation  is  incorrect. 


Tince  programming  difficulties  have  occurred  in  the  evaluati  t. 
of  "q.  (4j) , only  Eq.  (42)  will  be  compared  with  the  experimental 
data.  Drawings  AS-70-HT4  and  AS-70-1175  are  two  comparisons  of  theory 
with  experiment.  There  are  two  ways  these  predictions  can  be  improved: 
(l)  by  using  the  full  series  form  of  Eq.  (4y)  instead  of  just  the 
asymptotic  form  given  by  Eq.  (42),  and  (2)  by  a more  accurate  evalua- 
tion of  the  slope  correction  term  F.  The  expression  for  F given  in 
this  report  is  only  the  evaluation  of  the  true  slope  correction  term 
(-sin0o  + i^cos 6 ) at  the  stationary  phase  point  ' = -V^/V^.  The 

resulting  form  for  F is  then  taken  out  of  the  scattering  integral.  A 
much  more  accurate  way  to  predict  backscattered  results  would  be  to 
keep  everything  under  the  integral.  When  the  above  corrections  are 
included  in  the  form  for  the  backscatter  coefficient,  tb.  comparison 
of  theory  with  experiment  should  be  improved. 


Fresnel  high  Frequency  Backscatter 

The  intensity  equation  from  which  the  high  frequency  backscatter 
predictions  will  be  made  was  derived  in  the  Final  F.eport  under 
•'ontraet  1J00024-69-C-1  7‘>  and  has  the  following  form: 
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SURFACE  I 
FREQUENCY  = 500  kH 
h = 0.35  x 0.87 


BACKSCATTERING  COEFFICIENT  AS  A FUNCTION  OF  GRAZING  ANGLE 


If  the  amplitude  factor  is  disregarded,  then  Pq  is  simply  the 
directivity  pattern  of  the  source.  In  the  final  report  mentioned 
previously  and  in  the  previous  progress  report  under  this  contract, 
the  necessity  of  using  a directivity  pattern  and  the  Fresnel  phase 
approximation  was  discussed.  Also,  an  analytically  tractable  approxi- 
mation for  the  directivity  pattern  was  presented.  The  detailed 
results  for  the  intensity  calculation  based  on  the  assumption  of 
Gaussian  surface  statistics  were  presented  in  the  Final  Report  under 
Contract  N00024-69-C-1275.  However,  at  that  time  only  the  forward 
scatter  results  were  presented.  It  was  clear  then  that  Gaussian 
statistics  were  inadequate  for  backscatter  prediction  since  they  do 
not  predict  the  Lambert's  Law  type  behavior,  which  occurs  at  small 
grazing  angles;  consequently,  an  investigation  of  other  types  of 
probability  distributions  of  heights  was  undertaken. 


Consider  first  how  the  stochastic  nature  of  the  problem  is  taken 
into  account.  The  principal  quantity  of  interest  in  Eq.  (U5)  is  the 
bracketed  expression 
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If  the  processes  £ and  £ ' are  assumed  to  be  homogeneous  and  the 
following  conventional  approximations  are  made, 


e • e =0,  e • e_  = cos6  , e -e.  = sin6  , 7 = sinS  + sine, 
x 1 yl  r z 1 r r i 

(>+7) 
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then  according  to  Appendix  B of  the  Final  T-eport  under  Contract 
:000'4-69-  Eq.  (46)  has  the  equivalent  expression 


" Ein  6r\  7 + 


cos  6 -s2 


?^<> 


where  /^is  a natural  notation  for  . In  actually 

performing  the  derivatives  indicated  in  Eq.  (48),  it  will  he  assumed 
that  the  processes  > and  £’  are  also  stationary;  this  implies 
8h/8x.  = 0,  where  h is  the  rms  value  of  the  surface  height.  In 

l 

addition,  for  a stationary  process  the  correlation  function  depends 
only  on  the  difference  of  the  coordinates  i x-x’  and  q = y-y' . 
resulting  in  a considerably  easier  integration. 

The  most  interesting  feature  of  Eq.  (48)  is  that  only  knowledge 

of  the  characteristic  function  and  the  correlation  function  is  needed 

to  determine  the  quantities  in  Eq.  (46)  completely.  That  is,  the 

probability  distribution  function  that  is  required  is  a>( . and 

not  a>( t , f ' , C r ", ' ) ; as  a result,  the  Droblem  is  much  simpler, 

x x’  y 3y' 

The  correlation  functions  that  will  be  of  interest  in  this  prosy 
report  .are  the  following: 


hAUSCTAN  : C(  ? ,rj ) e 


E:  i'OEEETIAL:  C ( S , tj  ) e 


HAWKE! : C(F,  ,q) 


L +f.  +q 


It  was  pointed  out  in  Section  B that  the  exponential  correlation 
suffers  from  several  disadvantages,  most  notably  a lack  of  analiticity. 
In  view  of  this  problem,  a search  was  made  for  a function  which  would 
behave  like  an  exponential  function  moderately  near  the  origin,  and 
yet  still  be  an  analytic  function;  the  resulting  function  we  have 
named  the  Hawker  correlation  function.  Therefore,  in  the  actual 
calculation  of  backscattering,  only  the  Gaussian  and  the  Hawker 
correlation  functions  will  be  used. 

It  has  previously  been  shown  that  the  high  frequency  limits 
(or  asymptotic  forms)  of  the  intensity  are  obtained  by  expanding  the 
correlation  function  in  series  and  retaining  only  the  first  two  terms. 
Thus,  for  the  Gaussian  and  Hawker  correlation  functions  the  high 
frequency  limits  are  obtained  respectively  from 


.2  2 

C(?,r])  = 1 - i-EL  + 

G L" 


C(?,t))h  = 1 - — ^ + 


So  that  no  confusion  results,  it  is  necessary  to  keep  in  mind 
that  the  parameter  L,  which  is  the  correlation  length,  takes  on 
different  values  in  the  Gaussian  and  Hawker  correlation  functions. 


In  addition  to  the  two  correlation  functions,  three  different 
bivariate  statistics  will  also  be  considered:  Gaussian,  Exponential, 

and  Three-halves  (Student's  t with  degrees  of  freedom).  The  details 
of  the  characteristic  function  calculations  were  presented  in 
Appendix  A of  the  Final  Feport  under  N000P4-69-C-1  . so  they  will  be 

dispensed  with  here.  In  fact,  most  analytical  details  in  the  following 
calculations  will  be  omitted  so  that  only  the  final  results  will  appear. 


'aucsian  Bivariate  Distribution 


f; 


The  characteristic  function  for  the  Gaussian  bivari ate 


li  -tribution  is 


A-ik7(C-?’)^  = e-g(i-c) 

"pon  substitution  of  Eq.  (52)  in  (54),  one  obtains 


( 5) 


where  g = 4 sin  0 i V . W hen  Eq.  (55)  is  used  in  conjunction  with 
Eq.  ( -fd) . the  Fresnel  phase  approximation,  and  beam  function  approxi- 
mation in  Eq.  (45),  the  following  result  is  obtained: 
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(56) 


•where  s is  the  rms  slope  and  is  given  by 

hVT 

s = ~ j? 

and 


l6:rtKr,  r 

lo  00 

Her*  A is  the  in sonified  elliptical  area,  and  Y is  a constant  relat*  t 
to  the  beam  function . This  result  will  "be  iiscuss*  i in  • r.junctioi  wit! 
th<  results  for  the  other  bivariate  li rtributiunr . 


in.  rivariate  Distribution 


i'he  characteristic  function  for  the  Exponential  bivariate 
distribution  is 


-ik7(C-C')' 
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uhich  becomes  upon  substitution  of  Eq.  (^9), 


(99) 


(60) 


Using  Eq.  (60).  the  backscattered  intensity  is  obtained  in  a manner 
identical  to  the  Gaussian  calculation. 
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where  s and  0 are  given  by  Eqs.  (57)  and  (;  J). 


The  high  frequency  version  of  the  characteristics  function 
for  the  Three-halves  bivariate  distribution  is 


(5%2) 


(62) 


and  from  this  the  backscattered  intensity  is  obtained  as 
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It  should  be  noted  that  the  calculations  leading  to  the 
intensities  given  by  Sqs.  (56),  (6l),  and  (65)  were  performed  with 
the  Gaussian  correlation  function.  Comparison  in  Eqs.  (52)  and  (5;) 
shows  that 


(64) 


Hence  all  of  the  results  for  the  intensity  calculated  for  the  Gaussian 

correlation  function  can  be  extended  to  the  Hawker  correlation 

function  by  using  the  simple  substitution  given  by  (64).  The 

values  of  L_  and  L that  give  the  best  theoretical  fits  to  +he 
G H 

experimentally  measured  correlation  function  are  i.4c  in.  and  1.14  in., 
respectively.  The  slopes  in  each  of  the  cares  then  become 
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Comparing  Eqc.  (64)  and  (65),  it  appears  that  in  the  high 
frequency  limit  a surface  characterized  by  a Hawker  correlation 
function  behaves  as  a rougher  surface  than  one  characterized  by  a 
laussian  correlation  function  because  of  the  greater  rms  slope  ) 

This  is  an  interesting  feature  since  in  the  low  frequency  limit,  the 
opposite  is  true  because  the  Hawker  function  has  much  greater  correla- 
tion than  the  Gaussian  function  in  the  tail. 

In  addition  to  the  functional  similarity  obtained  when  using 

Lther  a Hawker  or  a Gaussian  correlation  function,  there  are  several 

other  interesting  features  inherent  in  the  high  frequency  limit. 

first  notice  that  there  is  absolutely  no  frequency  dependence  in  the 

intensities  given  by  Hus . ( 36),  (6l) , and  (63).  Hex'  , notice  that 

there  is  no  height  dependence;  rather  the  roughness  dependence  is 

"based  strictly  on  the  rms  slopes.  Finally,  it  is  3lear  from  th< 

parameter  G that  the  intensity  is  area  dependent  an  i that  the  ran  ~u 

dependence  is  given  by  1 r,  r and  not  l\  geometrical  acoustic.  , i.  •. 

- lo  00 

1 (r  +r  ).  In  previous  reports . it  ha.  been  shewn  th  i for  : lanar 
lo  00 

surfaces  the  geometrical  acoustics  range  Impendence  1 ( r.  +r  ' is 
btained  as  predicted  by  image  theory,  and  th  r is  n ins  nified 

area  iependence  for  a Lnfinit*  plane . , in  th*  inal  R* inde 

:ontract  N *-6§M?-l  ’ was  noted  1 at  for  n 1 rat*  1 u : 
surfaces  at  medium  frequencies  the  ran,  e f-p>..-ndenc<  was  rnasl- 
;e  met  si  cal  acou*  tics  . and  there  was  Litl  area  lependei  . rhus, 
in  passing  to  the  higl  fr<  piency  limit,  the  Fraunl  fer  ran  ■ iependenc 
which  has  been  characteristic  of  all  jther  • h*  ri*  s,  is  ;1*  arl> 
obtained.  However,  there  must  obviously  V some  transition  region 
which  none  of  the  other  th*;  rie:  hav*  tah*  • ini  ac  • mnt  . It  has  be*  n 
claimed  in  some  papers  that  even  in  th*  '•  fr*  U<  n limit  ' ••••.,■•  . 

for  example)  the  Fraunhofer  range  dependent*  is  >M.aiv  i. 


( 


The  theoretical  and  experimental  quantities  compared  in 
Hwgs.  AS-70-1194,  AS-70-1195;  and  AS-7O-II96  are  the  scattering 
coefficients  defined  ty  Eq.  (36): 
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(67) 


where  the  theoretical  hackscattered  intensities  given  by  Eqs.  (56), 

(6l),  and  (63)  have  been  used.  Care  was  taken  in  the  selection  of 
the  experimental  data  to  be  sure  that  they  actually  corresponded  to 
the  high  frequency  limit.  It  is  fairly  certain  that  the  data  shown 
for  surfaces  III  and  I are  in  the  high  frequency  limit,  since  the 
200  kHz  and  ^00  kHz  lie  quite  close  together;  that  is,  they  exhibit 
no  frequency  dependence.  Another  interesting  aspect  of  the  drawings 
is  the  notable  lack  of  agreement  between  experiment  and  the  theoretical 
predictions.  This  was  unexpected  since  the  Three-halves  and  Exponential 
distributions  resulted  in  extremely  good  agreement  between  experiment 
and  theory  for  the  coherent  intensity  and  for  the  forward  scattered 
intensity.  One  of  the  most  prominent  aspects  of  the  Three-halves 
theory  is  the  upturn  in  the  predicted  backscattered  intensity  at 
grazing  angles  below  40  deg.  If  shadowing  theory  had  been  incorporated, 
this  feature  would  not  be  so  pronounced  even  though  it  would  still  be 
present.  The  reason  for  this  behavior  is  that  at  low  angles  the 
Three-halves  distribution  predicts  that  the  major  intensity  component 
will  be  given  by 
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This  type 
inclusion 


of  dependence  at  low  grazing  angles  is  due 
of  slope  correction  terms  in  the  intensity 


strictly  to 
integrals . 


the 
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Examination  of  the  stationary  phase  backscattered  intensity  expression 
[Eq.  (42)]  does  not  reveal  a similar  difficulty  since  the  slope  was 
evaluated  at  the  stationary  phase  point  in  Eq.  (42). 


Because  of  the  lack  of  agreement  in  the  theoretical  and 
experimental  backscattering  expressions,  additional  aspects  on  the 
backscattering  problem  will  have  to  be  examined.  In  the  next  section, 
for  example,  a qualitative  investigation  of  composite  surface  scattering 
will  be  given.  Also,  it  is  hoped  that  in  the  next  progress  report 
the  integrations  of  certain  other  promising  bivariate  distributions 
will  be  completed,  and  some  results  for  composite  surfaces  will  be 
available. 

D.  Composite  Surface  Theory 

Recently  there  have  been  several  attempts  to  explain  backscattering 
effects  by  using  the  composite  surface  model.  However,  as  pointed 
out  in  the  introduction,  these  attempts  have  not  given  physically 
consistent  predictions  for  both  backward  and  forward  scattering.  In 
this  section,  qualitative  arguments  will  be  developed  which  show 
that  the  composite  nature  of  the  surface  can  have  a significant  effect 
only  on  the  backscattered  portion  of  the  intensity.  The  small  scale 
roughness  does  not  significantly  affect  the  specular  or  the  forward 
scattered  intensity. 

Assume  the  random  surface  process  is  represented  by 

r.(x,y)  = ^(x.y)  + '„(x,  y)  , (69) 

where  r and  * are  independent  stationary  random  functions  with 
zero  mean,  each  with  its  own  statistical  distribution  and  correlation 


■'I 


function.  In  addition,  the  processes  are  such  that  the  rms  values  of 
the  heights  obey: 


h^  » hOJ  and  \ > g . (7°) 

Equation  (70)  merely  implies  that  £ represents  some  large  scale 
roughness  -while  represents  some  small  scale  roughness.  Here  \ is 
the  source  wavelength,  which  must  be  about  the  same  order  of  magnitude 
or  greater  than  h^. 


Recall  that  for  the  pressure,  the  statistical  quantity  of  interest 
is  given  by 
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(71) 


where  tj.  is  defined  as  dC^/dy  and  is  independent  of  Based  on  this 

independence,  Eq.  (77)  can  be  written  as 


siner()-  coser^L)  + (r]2))(  ) , 


(72) 
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■■.here  ^ ^ represents  {c 
astit 

in  the  same  manner,  then  essentially 


lir.  the  brackets  mean  the 

stochastic  average,  if  it  is  assumed  that  g and  are  distributed 


(73) 


where  X is  the  general  notation  for  the  characteristic  function. 


However,  from  Eq.  ('70)>  one  sees  that 


< ) - 42?2(hi+i4)]  * 42-/2i>tl  • <7“) 

Clearly,  the  small  scale  roughness  is  insignificant  in  the  character- 
istic function. 


Now  examine  the  averages  tj  and  T,,.  If  shadowing  is  ignored, 
then  since  £ and  were  zero  mean  random  processes,  it  is  known 
that  tj^  and  are  also  zero  mean  random  processes  if  £ ^ and  are 
homogeneous.  However,  some  very  crude  shadowing  approximations  can 
be  made.  First  recall  that  tj  is  the  tangent  to  the  surface  at  the 
point  (x,y)  on  the  surface.  Then  if  the  possibility  of  the  point  (x,y) 
on  the  surface  being  shadowed  by  some  other  near  lying  surface  point 
is  ignored,  it  can  be  claimed  that  the  point  is  not  shadowed  if 


- tan£h  < Tj  < tanfi^  . (75) 

Hence  the  averages  of  ^ and  tj,__  are  not  given  by 


but  rather  by 


tan0 

(\)  / 

-tan0 


(76) 


(77) 


since  only  those  values  of  tj  that  are  not  shadowed  with  respect  to 
both  the  source  and  receiver  can  contribute  to  the  pressure  integral . 
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In  Eqs.  (76)  and  (77),  °(ri^)  is  the  probability  distribution  of  Tj 

and  -will  be  required  to  be  an  even,  zero  mean  function.  With  this  in 

mind,  various  cases  of  Eq.  (77)  can  now  be  investigated.  In  the 

specular  direction  tan0^  = tanS^,  then  <tj^>  and  <i],>  equal  zero  from 

Eq.  (77)-  Also,  in  forward  scatter  tan0^  does  not  differ  appreciably 

from  tane^  for  most  cases,  and  again  <r)^>  and  <t)?>  are  essentially 

zero.  However,  for  backscattering  tanQ^  ->  00 , and,  depending  on  the 

value  of  tan0^,  the  averages  can  differ  appreciably  from  zero.  For 

normal  incidence,  tan0  and  again  <r)  > and  <r),  > are  zero,  but  as 

r 1 

the  grazing  angle  approaches  zero,  Eq.  (77)  becomes 


<v 
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(78) 


This  is  the  case  where  <q^>  and.  <t]0>  take  on  their  greatest  value. 
Suppose  that  while  the  rms  value  of  the  heights  of  the  large  scale 
roughness  satisfies  h^»h^  as  can  often  be  the  case,  the  rms  value  of 
the  slopes  of  the  small  scale  roughness,  s-,  satisfies  s , »s ^ . Then 
<q ,->  S <q^>,  and  if  h^»h,  , then  Eq.  (72)  becomes 


sin0  x[ki>2h(r 
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A good  physical  example  of  a composite  process  where  this  would  be 
true  would  be  sand  dunes  made  up  of  angular  grains  of  sand.  The 
heights  of  the  dunes  would  be  very  much  greater  than  the  grain  sizes, 
but  the  slopes  of  the  angular  faces  of  the  grains  can  be  very  much 
greater  than  the  slopes  of  the  dunes.  A realistic  example  of  scat- 
tering by  this  type  of  composite  surface  is  the  scattering  of  light 
by  seemingly  flat  snow  fields. 


7 4 


Examination  of  Eq.  (79)  reveals  the  very  interesting  behavior 
that  can  he  exhibited,  by  certain  types  of  composite  surfaces.  For 
most  forward  and  for  all  specular  scattering,  the  lead  term,  which 
depends  only  on  the  large  scale  roughness,  dominates.  Also,  for  back- 
scattering  at  near  normal  incidence,  the  lead  term  dominates;  that 
is,  the  large  scale  roughness  is  again  the  major  factor.  However,  for 
backscattering  at  low  grazing  angles  the  second  term  in  Eq.  (7b)  is 
predominant.  Hence,  it  is  only  at  low  grazing  angle  backscatter  that 
the  small  scale  roughness  can  contribute  and  only  then  if  it  has  an 
appreciable  rms  slope  value. 


Let  us  summarize  the  preceding  arguments.  For  a composite 
surface  model  to  be  useful  the  following  conditions  must  be  met: 


hl  >>  h2 


X > h. 


s„  » s 
e:  1 


The  second  condition  merely  implies  that  the  high  frequency  limit 
must  not  have  been  reached  since  the  scattered  intensity  is  then 
strictly  slope  dependent,  and  the  small  scale  roughness  will  dominate 
the  scattering  effects  at  all  angles.  Finally,  for  physically  con- 
sistent results  to  be  obtained  for  both  forward-  and  backscattering, 
shadowing  considerations  must  be  employed. 
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